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j Mass  fraction  of  species 
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Dj  Diffusion  coefficient  of  species  j 
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ji-h  species 
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Abstract 

The  effects  of  flow  property  fluctuations  on  the 
stagnation  point  heat  transfer  rate,  as  calculated  by  the 
Fay-Riddell  equation,  are  investigated.  Fluctuations  in 
free  stream  temperature  and  density,  wall  temperature  and 
density,  and  stagnation  pressure  are  considered.  The  flow 
parameter  perturbations  are  included  in  the  Fay-Riddell 
equation  by  approximating  the  heat  transfer  rate  by  a 
Taylor  Series  truncated  to  second  order  and  expanded  about 
the  mean  values  of  the  five  independent  variables.  Taking 
the  time  average  of  this  expression,  the  first  order  terms 
in  the  fluctuations  drop  out,  and  an  expression  for  the 
time  averaged  stagnation  point  heat  transfer  rate  is  obtained. 
The  second  derivatives  in  the  Taylor  Series  expansion  are 
obtained  analytically  and  checked  numerically.  The  effects 
of  flow  property  fluctuations  on  the  Fay-Riddell  equation 
are  investigated  in  general  for  a wide  range  of  steady 
flow  conditions  and  fluctuations.  A comparison  between  the 
laminar  steady  flow  heat  transfer  rate  and  the  heat  trans- 
fer rate  with  fluctuations  shows  that  neglecting  the  fluctua- 
tions in  the  five  independent  variables  can  lead  to  errors 
in  the  calculated  heat  transfer  rate.  Of  the  five  variables 
considered,  the  fluctuations  of  the  temperature  at  the 
edge  of  the  boundary  layer  have  the  greatest  effect  on  the 
heat  transfer  rate.  Fluctuations  of  half  the  initial  value 
of  this  variable  result  in  a 5 to  10  percent  change  in  the 
heat  transfer  rate  over  the  range  of  initial  conditions 
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considered.  The  error  increases  significantly  with  the 
magnitude  of  the  fluctuations.  For  example,  a 100  percent 
fluctuation  in  this  temperature  results  in  a 15  to  35 
percent  error  in  the  heat  transfer  rate  for  the  range  of 
temperature  considered.  The  results  of  the  cases  investi- 
gated are  shown  in  graphical  form.  A computer  program 
implementing  the  solution  procedure  is  included. 


I . Introduction 


Background 

The  Re-entry  Nose  Tip  (RENT)  Test  Leg  of  the  5 0 Mega- 
watt Facility  of  the  Air  Force  Flight  Dynamics  Laboratory 
(AFFDL)  is  being  used  as  a ground  test  facility  for  evalua- 
ting the  performance  of  thermal  protection  systems.  It 
consists  of  a high  pressure,  high  voltage  arc  heater,  inter- 
changeable supersonic  nozzles,  a linear  motion  model  and 
probe  carriage,  an  exhaust  system,  and  a high  speed  Ambilog 
data  acquisition  system.  A schematic  of  the  major  system 
components  is  shown  in  Figures  1 and  2.  The  RENT  Test  Leg 
produces  a high  enthalpy,  high  pressure  airstream  ideally 
suited  to  the  study  of  ablation  phenomena.  Other  high 
pressure,  high  heating  rate  tests  dealing  with  transpira- 
tion cooling,  and  internally  water  cooled  leading  edges  can 
also  be  accomplished  (Ref  11:12). 

Calibration  of  the  RENT  Test  Leg  has  been  attempted  on 
several  occasions  by  personnel  from  AFFDL  as  well  as 
Aerotherm  Corporation.  During  the  calibration  attempts, 
null-point  calorimeters  were  swept  across  the  test  jet  to 
obtain  the  heat  flux  profiles  at  essentially  fixed  arc 
heater  conditions.  Large  fluctuations  in  heat  flux  have 
been  observed  for  all  nozzles  used  in  all  calibration 
attempts.  The  fluctuations  in  heat  flux  have  appeared  from 
run  to  run,  probe  to  probe,  as  well  as  for  one  probe  when 
swept  through  the  flow  from  a large  nozzle.  Fluctuations 
of  this  type  have  also  been  observed  in  other  arc  heated 
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RENT  Test  Leg  Schematic 


ARC  HEATER 


Model  Support 


facilities  used  for  nose  tip  testing  (Ref  2:31). 

There  exist  several  mechanisms  within  the  facility 
arc  heater  which  could  be  the  source  of  the  heat  transfer 
rate  variations.  Some  possibilities  are  the  fluctuations 
in  arc  heater  supply  voltage,  fluctuations  in  arc  heater 
supply  current,  variations  in  the  high  pressure  air  supply, 
rotational  frequency  of  the  arc  attachment  points,  and  the 
arc  shunting  frenquency.  More  extensive  studies  are 
required  before  a conclusion  can  be  made  as  to  the  definite 
source  of  the  fluctuations. 

Purpose  of  Study 

Although  the  exact  source  of  the  fluctuations  is  un- 
known, investigation  into  the  effects  of  flow  property  pertur- 
bations on  the  laminar,  steady  state,  Fay-Riddell  equation, 
which  is  presently  used  to  calculate  heat  transfer  rates, 
is  of  current  interest.  During  calibration  and  ablation 
tests,  the  AFFDL  personnel  measure  the  stagnation  point 
heat  transfer  rate  using  null-point  calorimeters  with  the 
slug  located  at  the  stagnation  point.  Also,  during  both 
calibration  and  test  runs,  the  stagnation  pressure  is 
measured  and  appears  to  have  significant  fluctuations. 

An  average  centerline  stagnation  pressure  and  an  average 
centerline  heat  transfer  rate  are  then  computed. 

Interpretation  of  ablation  nose  tip  testing  requires 
a knowledge  of  the  total  enthalpy  of  the  flow.  Currently, 
the  theory  of  Fay  and  Riddell  is  being  used  to  relate  the 
measured  heat  flux  and  stagnation  pressure  to  the  total 


enthalpy.  The  Fay-Riddell  theory  is  also  used  to  evaluate 
other  flow  parameters  such  as  temperature  and  density. 

The  validity  of  applying  the  laminar,  steady  state, 
Fay-Riddell  equation  to  a facility  which  has  significant 
fluctuations  in  flow  parameters  such  as  the  RENT,  is  a 
question  of  great  concern  to  the  AFFDL  engineers.  In  this 
study,  the  effects  of  fluctuations  of  various  flow 
parameters  on  the  stagnation  point  heat  transfer  rate  are 
investigated  in  an  attempt  to  determine  the  validity  of 
applying  the  Fay-Riddell  equation  to  a flow  field  where 
fluctuations  exist. 

Scope 

The  effects  of  flow  parameter  fluctuations  on  the  Fay- 
Riddell  equation  are  investigated  in  general  for  a wide 
range  of  flow  conditions  and  fluctuations  that  could  apply 
to  a large  number  of  test  facilities.  Then  a more  in 
depth  investigation  is  accomplished  using  specific  flow 
conditions  and  realistic  ranges  of  fluctuations  that 
actually  exist  in  the  RENT  Test  Leg. 

The  effects  of  fluctuations  in  free  stream  temperature 
and  density,  wall  temperature  and  density,  and  stagnation 
pressure  are  considered.  These  effects  are  investigated 
for  two  limiting  cases.  In  the  first  case,  the  fluctuating 
components  are  assumed  to  be  uncorrelated,  and  in  the  second 
case,  total  correlation  is  assumed.  For  purposes  of  this 
investigation,  the  heat  transfer  rate  is  assumed  to  be 
insensitive  to  the  boundary  layer  chemical  reactions.  The 
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air  is  assumed  to  be  composed  of  oxygen  molecules  and 
nitrogen  molecules.  However,  the  study  is  performed  in  a 
general  manner  to  allow  for  inclusion  of  dissociation 
effects  in  a later  study. 

Outline  of  Report 

This  report  is  written  to  describe  the  mathematical 
modelling  procedure  used  to  determine  the  effects  of  flow 
property  fluctuations  on  the  Fay-Riddell  equation.  The 
solution  procedure  described  in  this  report  is  implemented 
in  a computer  program  listed  in  Appendix  B. 

A brief  discussion  of  the  Fay-Riddell  theory  and  how 
it  is  used  to  determine  the  heat  transfer  rate  is  given  in 
Chapter  II.  Additionally,  a discussion  of  the  Taylor 
Series  expansion  method  used  to  mathematically  include  the 
flow  property  fluctuations  is  given  in  this  chapter. 

Chapter  III  provides  the  actual  mechanics  involved  in 
developing  the  mathematical  model  of  the  Fay-Riddell 
equation.  The  expressions  used  to  calculate  the  parameters 
in  the  equation,  such  as  the  transport  properties,  are  also 
included.  Additionally,  a description  of  the  steps 
involved  in  developing  the  Taylor  Series  expansion  and 
the  coefficients  for  this  series  are  shown  in  this  chapter. 

Chapter  IV  includes  the  results  of  the  cases  investi- 
gated to  determine  the  effects  of  flow  property  fluctuations 
on  the  Fay-Riddell  equation. 

Chapter  V contains  conclusions  and  recommendations 
resulting  from  this  study. 
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II . Theory 

Fay-Riddell  Theory 

A detailed  investigation  of  laminar  stagnation  point 
heating  with  dissociation  effects  was  reported  by  J . A. 

Fay  and  F.  R.  Riddell  in  1958  (Ref  4:78-85).  They  carried 
out  a numerical  solution  of  the  differential  equations  of 
the  laminar  stagnation  point  boundary  layer  using  tabulated 
values  of  the  properties  of  high  temperature  air  in  dissocia- 
tion equilibrium.  The  results  of  this  study  led  to  an 
equation,  referred  to  as  the  Fay-Riddell  equation,  that  is 
widely  used  to  calculate  the  stagnation  point  heat  trans- 
fer rate  of  blunt  bodies  in  hypersonic  flow. 

The  boundary  layer  equations  were  developed  in  general 
for  high  speed  flight  where  the  external  flow  was  in  a 
dissociated  state.  The  effects  of  diffusion  and  of  atom 
ionization  in  the  boundary  layer  were  included. 

For  purposes  of  their  analysis,  certain  assumptions 
were  made  to  establish  the  flow  condition  about  the  body. 

The  gas  was  assumed  to  be  a mixture  of  thermally  perfect 
gases  (i.e.  the  gas,  as  well  as  each  component,  obeys  the 
ideal  gas  law).  The  flow  in  the  boundary  layer  was 
considered  to  be  laminar  and  the  boundary  layer  thickness 
small  compared  to  the  curvature  of  the  body.  In  addition, 
mass  diffusion  due  to  pressure  gradients  was  neglected. 

The  local  heat  transfer  rate  to  the  body  was 
determined  by  the  sum  of  the  conductive  and  diffusive 
transports.  Thus 
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where  the  x and  y directions  are  tangential  and  normal  to 
The  body,  respectively.  Fay  and  Riddell  then  introduced 
dimensionless  temperature  and  enthalpy  distributions  and 
simplified  the  previous  expression  to  obtain  an  equation 
for  the  stagnation  point  heat  transfer  rate.  It  was 
expressed  as 


V 


where 


(2) 


Nu=  % Y.  Cpw  Re  = lie  x (3) 

Kw(h3-hJ  ; V* 

Then  non-dimensional  forms  of  the  continuity,  momentum, 
and  energy  boundary  layer  equations  were  solved  to  derive 
a simplified  expression  for  Nu/VRe  for  both  the  frozen  and 
equilibrium  boundary  layers.  This  investigation  deals  only 
with  the  results  of  the  equilibrium  boundary  layer.  First, 
an  expression  was  obtained  for  the  Nu/\ARe  ratio  assuming  a 
Lewis  Number  of  unity.  It  was 


JSIu  _ Q.67/W-*  (4) 

\/R&  Ww 

Then  for  other  values  of  the  Lewis  Number,  a correction  of 
the  form 


(n“//rI)u,i  ihj 


± 


was  included  where 


h„  = Z cis  ( - h° ) (6) 

Combining  these  two  expressions,  an  equation  was  obtained 
for  the  ratio  of  the  Nusselt  Number  over  the  square  root  of 
the  Reynolds  Number. 

/ \olf 

N u __  0.  (o~7 1 Ps  Us  1 
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Substituting  this  expression  into  the  equation  for  the 
stagnation  point  heat  transfer  rate 
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Fay  and  Riddell  then  assumed  Pr=0. and  simplified  the 
expression  to 


A recommendation  was  made  to  replace  the  factor  0.9^  with 
0.76  Pr  for  other  Prandtl  Numbers.  Incorporating  that 
factor  into  the  above  expression  resulted  in  the  followings 


^ - O 76 


(10) 


For  this  study,  a modified  Newtonian  flow  was  assumed; 


therefore,  the  stagnation  point  velocity  gradient  is 


for  a spherical  nose  in  hypersonic  flow.  The  free  stream 
pressure  is  much  less  than  the  pressure  at  the  edge  of  the 
boundary  layer  and vas  neglected.  Substituting  the 
expression  for  the  velocity  gradient  into  the  heat  trans- 
fer rate  equation  and  rearranging,  the  final  form  of  the 
Fay-Riddell  equation  used  for  this  analysis vas  obtained. 


Throughout  this  report,  q\fR  will  be  referred  to  by  Q. 

Fluctuation  Modelling  Theory 

A method  was  chosen  that  would  enable  incorporation  of 
the  flow  property  variable  fluctuations  into  the  Fay-Riddell 
equation.  The  free  stream  temperature  and  density,  the  wall 
temperature  and  density,  and  the  stagnation  pressure  were 


chosen  as  the  independent  variables.  For  purposes  of  this 

investigation,  the  five  independent  variables  are  expressed 

in  terms  of  a time  averaged  and  a fluctuating  component. 

Ts  = Ts  + Ts' 

(13) 

T*  - Tw  + Tv/ 

(14) 

?s  = Ja  * 9 a' 

(15) 

~ 9 w 

(16) 

p.  = p~  + p / 
r*z  Ktz 

(1?) 

The  barred  quantities  denote  average  values  and  the  primed 
quantities  denote  the  fluctuating  components. 
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Following  the  method  introduced  by  Weeks  (Ref  12),  and 

used  by  Lazdinis  (Ref  9)»  stagnation  point  heat  transfer 

rate  can  be  approximated  by  a Taylor  Series  truncated  to 

second  order,  and  expanded  about  the  mean  values  of  the 

five  parameters  T , T , p , p , and  ^t0.  The  resulting 

s w vs  2 

expression  for  the  instantaneous  value  of  the  heat  trans- 
fer rate  is  thus  a function  of  both  the  average  and  the 
fluctuating  terms  of  the  five  independent  variables. 

Q(Qs  ) Ft*  ) — pw;  ~Ts  ^ta)  4-  ps  4-  pw 

^Ps  ^Pw 

4-  j'  + ^0_  T/  + io.  p ' _L  second  order  terms  (18) 
^Ts  + 

The  derivatives  in  this  expression  are  evaluated  using  the 
average  values  of  the  independent  variables.  Taking  the 
time  average  of  this  expression,  the  first  order  terms  in 
the  fluctuations  drop  out,  and  an  expression  for  the  time 
averaged  stagnation  point  heat  transfer  rate  is  obtained. 
The  resulting  equation  used  for  this  investigation  is 


Q=  Q.(ft,P.,T„T,,pt.)+  + + Ij£  <Ty> 


, d*Q  /r'A  I d*Q  /p'2\.2  A2Q_  /'r'\  I 2 ,3*Q 


> 


,2^Q  / 'p  /v  , 2 /D  't'\  , 2 i!a  /PJj  ' 

f^Tw<esvv>  ^es^stl>  ^V?"r*'  + 

+ 2i Op.?''  P"' > + 2^<t:t‘>  + 2S&<TO + 2^.<t-'p'->]  ‘ 19  ’ 
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Ill . Mathematical  Formulation  of  the  Problem 

Thermodynamic  and  Transport  Properties 

The  mechanics  of  formulating  the  problem  involve 
applying  the  fluctuation  modelling  theory  to  the  form  of 
the  Fay-Riddell  equation  shown  in  the  previous  chapter. 

This  procedure  involves  writing  the  transport  properties 
and  parameters  in  the  equation  in  terms  of  the  five  indepen- 
dent variables.  The  required  derivatives  are  then  deter- 
mined analytically  to  allow  for  an  in  depth  investigation 
into  which  parameters  in  the  equation  are  significantly 
affecting  the  magnitude  of  the  derivatives.  Due  to  the 
complexity  of  the  equation,  the  derivatives  are  also 
determined  numerically  to  ensure  the  accuracy  of  the 
analytical  derivatives. 

In  calculating  heat  transfer  rates  at  high  speeds, 
variations  of  the  transport  and  thermodynamic  properties  j 

must  be  accounted  for.  Curve  fit  Prandtl  Number  data, 

| 

used  by  AFFDL  personnel,  is  used  in  this  study  (Ref  10). 

The  curve  fit  data  expresses  the  Prandtl  Number  as  a 
function  of  the  ratio  of  the  wall  enthalpy  and  a reference 
enthalpy.  If  h is  less  than  0.005»  the  Prandtl  Number 
is  set  equal  to  0.??.  For  hw/h-0  up  to  2.0,  the  Prandtl 
Number  is  computed  using  the  least  square  curve  fits. 

For  0.005<hw/ho£  0.1: 

Pr=0 . 8245251-0 . 1212104xl02(hw/ho)+0 . 349395xl03(hw/ho)2 
-0.4212513xl04(hw/ho)3+2.404883xl05(hw/ho)4 
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(20) 


-0.530703lxl03(hw/ho)5 

For  0.1  f h . /h  < 2.0: 

Pr=0. 788199-0. 1287464 (hw/ho)+0.23864lxl0-1(h w/hQ)2 

+0.102987xl0-1(hw/ho)3  (21) 

The  variation  of  Pr  with  wall  temperature  is  shown  in  Figure 

3- 

The  Lewis  Number  is  the  ratio  of  the  mass  diffusivity 
to  the  thermal  diffusivity.  Like  the  Prandtl  Number,  the 
Lewis  Number  is  also  calculated  by  curve  fit  data  (Ref  10). 
It  is  a function  of  the  ratio  of  the  enthalpy  at  the  edge 
of  the  boundary  layer  and  a reference  enthalpy.  If  h /h 

o U 

is  less  than  or  equal  to  0.1,  the  Lewis  Number  is  set  equal 

to  1.4.  If  h /h  is  greater  than  0.1,  the  following  least 
s o 

square  curve  fit  is  used: 

Le=l.  43691-0. 549-19(hs/ho)-0.07l6(hs/ho)2 

+0.06833(h  /h  )3  (22) 

s 0 

The  variation  of  Lewis  Number  with  T is  shown  in  Figure  4. 

s 

The  ratio  h^/h  is  also  calculated  from  a least 
u s 

square  curve  fit  of  data  as  a function  of  the  enthalpy  at 
the  edge  of  the  boundary  layer  (Ref  1).  The  following 
least  square  curve  fit  is  used  for  this  ratio: 

hn/h  =-0. 124554+0. 22780xl0-3h  -0 . 2944l7xl0_7h  2 
US  s s 

+0.184600xl0-11h  3-0.432111xl0_l6h  4 (23) 

s s 
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For  the  temperature  ranges  considered  in  this  investi- 
gation, the  viscosity  is  essentially  dependent  only  on  the 
temperature.  To  account  for  the  viscosity  variations  in 
the  high  speed  boundary  layer,  the  Sutherland  Formula  is 
used. 


U - Tr +198.6  f T \ (24) 

T + 198-6  V TR  / 

The  viscosity  of  air  for  temperatures  in  degrees  Rankine  is 
given  by 


M - 2 . 2 7 x 108  T (25) 

T + 198.6 

where  /u  is  in  slugs/ft-sec  (Ref  5*^9) • 

The  enthalpy  was  calculated  in  a form  that  would  allow 
inclusion  of  all  the  species  of  air,  although  only  oxygen 
and  nitrogen  molecules  were  considered  in  this  study.  The 
enthalpy  of  the  mixture  is  determined  by  taking  the  sum  of 
the  product  of  the  species  mole  fractions  and  the  correspond- 
ing species  enthalpy  (Ref  10:1?).  The  species  enthalpy  was 
calculated  from 


h,-h°  _ 5«2flvi)  , (<Vl)e,;  S»j  e T (26) 

T ' 2 T(e"^- 1)  r ,t  e‘Vr 

"t  h 

where  the  subscript  j represents  the  j species  and  the 

"t  h 

subscript  2 represents  the  2X  electronic  level  of  that 
species.  This  expression  includes  the  translational, 
rotational,  vibrational,  and  electronic  contributions  of 
each  species.  The  use  of  this  equation  requires  the 
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specification  of  molecular  constants  Oj  , , h/  , , 

and  for  each  species.  The  values  of  these  quantities 
for  oxygen  and  nitrogen  molecules  are  given  in  Table  I. 

Table  I 

Thermodynamic  and  Molecular  Constants 


Constant 


N, 


Atoms  Per  Molecule 


2.0 


2.0 


Characteristic  Vibrational 
Temperature  (°K) 


3.35324x103 


2.23897x103 


Enthalpy  of  Formation 
( cal/mole ) 


0.0 


0.0 


Electronic  Levels 


4.0 


5.0 


Level  I 


Degeneracy 


1.0 


3.0 


Energy  (cal/mole) 


Level  II 


0.0 


0.0 


Degeneracy 


3.0 


2.0 


Energy  (cal/mole) 


1.43685x10- 


2.26370x10 


Level  III 


Degeneracy 


6.0 


1.0 


IT 


Energy  (cal/mole) 


1.70475x10- 


3.77250x10 


Level  IV 


Degeneracy 


1.0 


3.0 


Energy  (cal/mole) 


1.71540x105 


1 .03198x105 


Level  V 


Degeneracy 
Energy  (cal/mole) 


3.0 

1.42390x103 


These  expressions  for  Prandtl  Number,  Lewis  Number, 
viscosity,  and  enthalpy  are  substituted  into  the  Fay-Riddell 
equation  to  express  the  heat  transfer  rate  in  terms  of  the 
five  independent  variables. 

Analytical  Derivatives 

Before  taking  the  desired  derivatives  analytically, 
the  Fay-Riddell  equation  is  rearranged  to  group  together 
terms  that  are  functions  of  each  independent  variable. 

Start  wi.th  the  Fay-Riddell  equation  in  the  form  of 

(K-  (27) 

0 2^ 

Let  C=0.?6(2)  * D and  rearrange  to  obtain 


Q - 0.76  Pr  V/-JW 


4r-i)h 

l~>£. 


-Ob  01  CHI  OIS  0.25  r ot 

Q-CPr  /iw  ps  Pt2  J fXs 


(l  Tz~\)hc 


r (hs-hw) 


(28) 


Now  define  the  following  group  of  terms  that  are  a function 
of  Tg  as 


o + 

^ = Ms 


■d 


0.52 

e 


-Oh* 

hs 


(29) 


and  the  group  of  terms  that  are  a function  of  Tw  as 


(3  = Pr 


-ob  o.l 

A*  w 


(30) 


Also  define 


ar  = 


_ 0.1  o.l 5 _ 0 2 5 

C Pw  Ps  Ptz 


(3D 


Rewriting  the  Fay-Riddell  equation  to  include  <x,  £3  , and  Jf 


1? 


(32) 


Q=  x(?s,pwjPtz)  ■ H5-°<(Ts) - (3(TJ-h*-  «(TS)J 

Simplifying  further,  let 

A = hs  <»<■  (Ts)  (33) 

and 

B = h„  @(TW)  (3^) 

The  equation  can  then  be  written  as 


Q = *[gA  - * B]  (35) 

where  ^ is  a function  of  the  density  at  the  wall,  the 
density  at  the  edge  of  the  boundary  layer  and  the  stagnation 
pressure;  ^ and  B are  functions  of  only  the  wall  temperature; 
and  A and  o.  are  functions  of  only  the  temperature  at  the 
edge  of  the  boundary  layer.  Using  this  form  of  the  Fay- 
Riddell  equation,  the  second  derivatives  required  for  the 
Taylor  Series  expansion  can  be  taken. 

The  second  derivative  of  the  heat  transfer  rate  with 
respect  to  the  wall  temperature  is 


&ZQ  _ 
AT/ 


A A*g>  oc  A*B 


AT 


2 

W 


atw 


(36) 


The  first  and  second  derivatives  of  0 are  required  in  this 
equation  and  also  in  some  of  the  other  analytical  derivatives. 


Aft  _ A ( Pr  M*  ) 

A Tw  A”fT 

•0.6  -0<?  . "<•<>  0.1 
Aft  - O.  I Pr  /aw  _ 0. 6 Pr  A Pr 

A Tv  A Tw 


(37) 

(38) 
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and 


, _ -Ckfc 

<a  3 — 0.  i R* 

a>T~ 


09  /-v  "(-9/x  \2  -/.6  -a<? 

a 3-*yy  -0.0?  R'  ,uw  / ^ I 0. 12  Fr  41*  <)m  a R- 


^T, 


O.l  I 


_0.%Pr  yUw  /IS. 


\2  r^"'6  O.l 

_ 0.6  B>  /U„ 


w/ 

2r 


^Pr 

<^TW2 


ai*/  yr„ 

(39) 


Calculation  of  the  (3  derivatives  requires  the  first 
and  second  derivatives  of  the  viscosity  and  the  Prandtl 
Number  with  respect  to  the  wall  temperature.  These 
derivatives  are  shown  in  Appendix  A. 

Additionally,  -^~z  . as  well  as  other  derivatives 

O « vv 

taken  later,  requires  the  first  and  second  derivatives  of  B. 


aB  _ 

— ( hw  b) 

(40) 

£TW  V V 

aB  . 

h>„  + ( 

3 

(41) 

^Tw 

<^TW 

and 


a2B  . 

- hw 

1 2 

^ + f 

3 a2k, 

(42) 

aV  ■ 

aV 

aV 

In  addition  to  the  ^ derivatives,  this  expression  requires 
the  first  and  second  derivatives  of  the  enthalpy.  These 
derivatives  are  given  in  Appendix  A. 

The  second  derivative  of  the  heat  transfer  rate  with 
respect  to  the  temperature  at  the  edge  of  the  boundary 
layer  is 


ik  B 

aTs2 


(43) 
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The  first  and  second  derivatives  of  «.  are  required  in  this 
expression  and  in  other  analytical  derivatives. 


_ 


*TS 

<>T5 


a_  { mV 


+ ( L e 0 


o.  6 

O f ps 


I + (Le  - 0 h 


A 0 lf 

^ 4-  A*s 


Le“ 


^T5 


+ 0.52yu^/v)Le°‘+8  ^Le 
\hs/  ^Ts 


and 


(44) 


(45) 


x2  -16 
6 ■*  ~ - 0.2H-  /JL% 


/ 0-52  lV 

+ (Le  - I)  hc 


<)ms\  + 0.4/6  ^6Le°"/. ha_)  Akl  Ails 


L ^ , 


-I- 


(Le-  I)  ifo) 

_f_  0.4yUs 

V(Le5l|)ky 

>TS 

o>Ts 

l K /_ 

<)TS  ^ts 


6^ 


I P.U.  Q“|  -O'*0  M v/h0\  ^ Of/,  -Of e 2 

q_l.04/xs  Le  <)  Le  <)(  nJ  0.  52  yUs  / bD  Le.  ^ L< 

6Ts  ^Ts  ( hs  / ^Ts; 

/u*(L/“l)  ,S*(£f)  _0.249^^7/i^\Le  48 /AU_\2 

^Ts2  (hs/  \^Tsj 


+ 


(46) 


Calculations  of  the  cx.  derivatives  require  the  first  and 
second  derivatives  of  the  viscosity,  hD/hs  ratio,  and  the 
Lewis  Number  with  respect  to  the  temperature  at  the  edge  of 
the  boundary  layer.  These  derivatives  are  shown  in 
Appendix  A. 

The  first  and  second  derivatives  of  A with  respect  to 
the  temperature  are  also  required  in  -5-^. 


6 A _ 6 ( hs  cx ) 

6TS  - yr; 


(A?) 
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A ^ ck 

<^Ts  <^)TS 


^ A _ 1^  ^ <*  _^_  2 <^s  ^ q<  , c*.  i^  h 

zV  ' ZV  ^Ts  ^Ts  aTs 


^Ts  IV  ^Ts  <^TS  ^Ts2 

The  second  derivative  of  the  heat  transfer  rate  with 
respect  to  the  wall  density  is 


a2<?  _ p-A  - o<-B 
<*e*2  ” <^e«2  L 


The  first  and  second  derivatives  of  X with  respect  to  the 
wall  density  are  needed  in  this  expression  and  other  cross 
term  derivatives. 


d*  _ _J_ 
d P»/  d Pw 

Y~  -O.C 


^>2*  _ 


0.2  5 o.l  o.  IS  0.2  5 


Ps  Pt4 

(5D 

0 25  -0  9 

Z 

(52) 

0 25  -1.9 

' t2  (V 

(53) 

The  second  derivative  of  the  heat  transfer  rate  with 
respect  to  the  density  at  the  edge  of  the  boundary  layer  is 


\ - ^L-  BA  - «■  B 

's  d Ps  l 


L*  J 

The  first  and  second  derivatives  of  X with  respect  to  the 
density  at  the  edge  of  the  boundary  layer  are  needed  in  this 
expression  and  other  cross  term  derivatives. 


\ v X r 0 25  0.1  O 15  0. 25 

Hr  = ik  P*,(2)  ^ p“  J 
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(56) 


~ = 0.13  6105  pw° ' pti° 25  es°  8S 
^ 2 - 0.  1 1 5 68 9 pj>'1  Ptz°  25  Ps"1’85  (57) 

<3CS 

The  second  derivative  of  the  heat  transfer  rate  with 
respect  to  the  stagnation  pressure  is 


£3  A - 01  B 


i 

(58) 


The  first  and  second  derivatives  of  X with  respect  to  the  ' 

1 

stagnation  pressure  are  needed  in  this  expression  and  other 
cross  term  derivatives t 

i 


b V 

= * \o.K>(2)ot%:\:'s  fC5l 

(59) 

6 * 

- 0.22&8lf\  e*'  P3°  ‘5  Pt;°75 

(60) 

ft 

- -0. 170(309  Vs  Pt^75 

(61) 

The  other  analytical  derivatives  required  in  the  Taylor 
Series  expansion  are  the  derivatives  involving  the  cross 
terms  of  the  five  independent  variables. 


32Q  _ 

t 

<)A  <)  £ 

()c* 

Ts  ^ 

<^TS  <^"0_ 

d*Q 

A - 

<* 

^ P*V 

L 

^T»v_ 

<^Q  - 

6 1 

~A  ^ - 

c)ps 

_ «)TW 

(62) 

(63) 

(64) 
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L 


< 


A Q _ ^ X A 1^ 

<^Pt2  _ 


<3* q _ b % ^A 


<^TS  & Rw  L*T. 


P - TV-  B 


d Q — ^ y A A D b ot,  n 

A~T^  Afs  «)  ps  _ATS  ATs 

— 

^ Q — ^ y c)A  q <^ <*  q 

<STs^f?z  <^Pt2  — ATs 

-AA_  - fe.  a - «*  ft) 

<^Pw<)Ps  ^ fs 


= O.OI36IO  ps° 8 5 Ptf 2 5pw0'  Y$  A - ■*  B) 


^ Q 2 — AL  /r -A- c*B) 

deUR*. 


: 0.022  68  tf  £XV’fp A - 01  B) 


^2Q-  ~ -i_  i*-(RA-<*B) 

APs^^t*  A ft  2 ^ Ps 


= 0.0340262  e°J  Pt°75pl°-9Sl 3A-*B)  (74) 


All  the  terms  on  the  right  side  of  these  equations  have 
been  previously  shown  in  this  chapter. 

The  second  derivatives  of  the  heat  transfer  rate 
mentioned  previously  were  also  checked  numerically  to 
ensure  accuracy  of  the  analytical  derivatives.  This  was 
accomplished  by  using  central  difference  formulas  (Ref  7:21- 
22).  A difference  of  less  than  one  percent  was  noted  between 

the  magnitude  of  the  analytical  and  numerical  derivatives 
for  the  second  derivatives  used  in  the  Taylor  Series 
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expansion . 


Development  and  Structure  of  the  Mathematical  Model 

The  model  is  developed  so  that  for  different  sets  of 
initial  conditions,  the  magnitude  of  the  fluctuations  of 
the  independent  variables  can  be  varied  and  all  possible 
combinations  of  independent  variable  fluctuations  can  be 
investigated.  The  initial  conditions  are  the  assigned 
average  values  of  the  five  independent  variables.  The 
effects  of  the  fluctuations  are  determined  by  comparing  the 
heat  transfer  rate  computed  using  only  the  initial  condi- 
tions to  the  heat  transfer  rate  computed  using  the  Taylor 
Series  expansion  shown  in  Equation  19. 

The  fluctuation  model  is  incorporated  into  a Fortran 
computer  program.  A listing  of  the  program,  description  of 
the  subroutines,  and  a listing  of  the  Fortran  variables  are 
given  in  Appendix  B. 

The  thermodynamic  and  transport  properties  are  first 
calculated  using  the  initial  values  of  the  five  independent 
variables.  Knowing  these  properties  for  a given  set  of 
initial  conditions,  the  heat  transfer  rate  without  fluctua- 
tions can  be  computed  using  the  form  of  the  Fay-Riddell 
equation  shown  in  Equation  12.  This  value  is  used  for 
comparison  purposes  after  the  heat  transfer  rate  with 
fluctuations  is  computed.  Also  from  the  initial  conditions 
the  values  of  <x  , (3  , 2 , A,  and  B can  be  calculated.  The 
expressions  given  in  Appendix  A are  then  used  to  calculate 
the  magnitude  of  the  first  and  second  derivatives  of  the 
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transport  and  thermodynamic  properties.  The  first  and 
second  derivatives  of  the  heat  transfer  rate  are  then 
computed  for  the  given  set  of  initial  conditions. 

The  fluctuation  components  of  the  five  independent 
variables  are  normalized  with  respect  to  the  initial 
values  of  the  respective  variables.  The  normalized  fluctua- 
tions are  represented  by 


£,= 

Ps' 

Ps 

(75)  £ = 

Tv/ 

Tw 

(78) 

P/ 

Pw 

(76)  £ = 

'■'5 

P*/ 

Ptz 

(79) 

V 

Ts 

(77) 

and  the 

cross 

term  derivatives  are  normalized  in 

a similar 

manner.  In  the  case  of  perfect  positive  correlation,  the 
following  relationship  between  all  fluctuating  variables  is 
assumed 

<£i  £,>  = (<  (80) 

Since  the  fluctuation  components  are  normalized,  the  second 
derivatives  in  the  Taylor  Series  expansion  must  be  multi- 
plied by  the  normalizing  quantity  to  allow  for  a consistent 
comparison  with  the  heat  transfer  rate  without  fluctuations. 

The  heat  transfer  rate  with  fluctuations  can  now  be 
calculated  by  including  the  second  derivatives  and  the  de- 
sired magnitude  of  the  fluctuations  considered.  The  percent 
difference  between  the  two  heat  transfer  rates  can  be 
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plotted  against  the  magnitude  of  the  fluctuations  of  the 
independent  variables  to  show  what  effect  different  combina- 
tions of  fluctuations  have  on  the  heat  transfer  rate  as 
calculated  by  the  Fay-Riddell  equation. 

A check  was  also  performed  for  each  case  to  ensure 
that  the  heat  transfer  rate  approximated  by  the  Taylor 
Series  expansion  was  accurate.  This  was  done  by  computing 
the  heat  transfer  rate  with  the  fluctuations  included  in 
the  value  of  the  independent  variables,  and  comparing  this 
result  to  the  full  Taylor  Series  expansion  including  the 
first  derivatives. 


> fw  ♦ pw  j '^st  Ts  |TW  j f t2  + j = Q0  ( Pn  ) 


+ n ' , £0.  D ' + ^ Q_  -r  ' + 7-  ' , iO_  p ' SfcoMQ  OR0£  R 

JTS'S  ZTW  rER*NS 


(81) 
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IV.  Results 


The  effects  of  flow  parameter  fluctuations  on  the 
Fay-Riddell  equation  were  investigated  in  general  for  a 
wide  range  of  steady  flow  conditions  and  fluctuations. 

The  following  set  of  initial  conditions  was  chosen  as  a 
baseline  for  all  cases  considered: 

ps  =0.005199  slugs/ft3 
p^=0. 003622  slugs/ft3 
Ts=7908.0  °R 
Tw=1170.0  °R 

Pt2=7597.°  lb/ft2 

These  are  typical  values  of  the  five  independent  variables 
obtained  from  RENT  Facility  data. 

Case  I 

The  first  case  investigated  involved  varying  each  of 
the  five  flow  parameters  separately  over  a wide  range  of 
initial  values  while  keeping  the  other  four  parameters  set 
to  the  baseline  values.  For  each  of  these  variations,  the 
changing  independent  variable  was  also  fluctuated  over  a 
range  of  S.  =0.2,  0.4,  0.6,  0.8,  1.0.  The  range  of  varia- 
tions of  the  initial  conditions  for  the  five  independent 
variables  was 

1x10"^  Ps  < lxlO-2  (slugs/ft3) 
lxlO~3<  < lxlO”1  (slugs/ft3) 

4000  sT  i 10,000  (°R) 

b 
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1000  5 Tw  f 4000  (°R) 

20000  -Pt2  - 100,000  (lb/ft2) 

The  results  of  these  cases  are  shown  in  Figures  5 through 
9 with  -the  magnitude  of  the  fluctuations  plotted  versus 
the  percent  error  in  the  heat  transfer  rate.  The  percent 
error  in  the  heat  transfer  rate  is  the  percent  difference 
between  the  heat  transfer  rate  without  fluctuations  (Q  ) 
and  the  heat  transfer  rate  with  fluctuations.  The  value  of 
Qq  varies  each  time  the  initial  conditions  are  changed  so 
it  is  given  in  the  figures  for  each  set  of  initial 
conditions  investigated.  Additionally,  the  parabolic 
equations  of  the  lines  are  shown  on  each  figure,  as  well 
as  the  coefficients  determined  for  each  set  of  initial 
conditions . 

Of  the  five  variables  fluctuated,  the  variation  of 
the  temperature  at  the  edge  of  the  boundary  layer  had 
the  greatest  effect  on  the  heat  transfer  rate.  The  percent 
error  in  heat  transfer  rate  increased  positively  as  the 
temperature  was  decreased.  For  £=0.6,  the  error  in  the 
heat  transfer  rate  was  5*3  percent  of  Qq  for  TS=10,000°F. 
and  increased  to  12.5  percent  for  Ts=4000°R. 

The  variations  of  the  density  at  the  wall,  density  at 
the  edge  of  the  boundary  layer,  and  stagnation  pressure 
gave  similar  results.  As  the  initial  values  of  those 
variables  were  increased,  the  percent  error  in  the  heat 
transfer  rate  decreased  in  magnitude.  However,  the  change 
in  the  percent  error  was  negligible  compared  to  the  change 
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Fig.  5-  Percent  Error  in  the  Heat  Transfer  Rate 
for  Various  Initial  Conditions  and  Fluctuations  of  Ts 
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Fig.  6.  Percent  Error  in  the  Heat  Transfer  Rate  for 
Various  Initial  Conditions  and  Fluctuations  of  p s 
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Fig.  8.  Percent  Error  in  the  Heat  Transfer  Rate  for 
Various  Initial  Conditions  and  Fluctuations  of  rt? 
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Fig.  9.  Percent  Error  in  the  Heat  Transfer  Rate  for 
Various  Initial  Conditions  and  Fluctuations  in  Tw 
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in  the  percent  error  for  variations  in  the  temperature 

initial  conditions.  For  £.=0.6,  the  error  in  the  heat 

transfer  rate  was  -2.3  percent  of  Q for  n =1.0xl0-^ 

O \ s 

slugs/ft-^  and  decreased  to  -2.2  percent  for  ^s=1.0xl0-2 
slugs/ft-^.  The  effect  of  wall  density  variations  was 
slightly  less  than  the  effect  of  variations  of  the  density 
at  the  edge  of  the  boundary  layer.  Also,  variation  of  the 
initial  stagnation  pressure  had  very  little  effect  on  the 
magnitude  of  the  error  in  the  heat  transfer  rate.  For 
£=0.6,  the  error  in  the  heat  transfer  rate  was  -3.4  percent 
for  pt2=20,000  lb/ft2  and  decreased  to  -3.2  percent  for 
Pt2=100,000  lb/ft2.  Figures  6,  7,  and  8 show  only  the 
upper  and  lower  limits  of  the  ranges  for  initial  conditions 
considered.  Other  values  within  each  range  were  investi- 
gated and  all  curves  lie  between  the  two  curves  shown  on 
each  graph. 

As  seen  from  Figure  9»  the  initial  value  of  the  wall 
temperature  had  an  effect  on  whether  fluctuations  of  this 
variable  would  increase  or  decrease  the  heat  transfer  rate 
and  to  what  extent.  The  variation  of  the  wall  temperature 
initial  condition  did  not  influence  the  heat  transfer 
rate  monotonically  as  was  shown  for  the  other  four 
independent  variables.  As  the  wall  temperature  baseline 
was  varied  from  1000°R  to  2000°R,  the  magnitude  of  the 
heat  transfer  rate  error  was  decreased.  Throughout  this 
range,  the  error  tended  to  decrease  the  heat  transfer  rate. 

For  T equal  to  2500°R,  the  magnitude  of  the  error  is  only 

w 
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2.3  percent  of  Qq  but  now  increases  the  heat  transfer  rate. 
For  Tw  equal  to  3000°R,  the  magnitude  of  the  error  is  very 
small  but  again  lessens  the  heat  transfer  rate.  As  the  wall 
temperature  baseline  was  varied  from  3000°R  to  4000°R,  the 
magnitude  of  the  heat  transfer  rate  error  was  increased  and 
tended  to  decrease  the  heat  transfer  rate.  The  second 
derivative  of  the  heat  transfer  rate  with  respect  to  the 
wall  temperature  has  been  shown  to  be  a function  of  the 
first  and  second  derivatives  of  ^ and  B.  The  erratic 
behavior  of  between  2000°R  and  4000°R  can  be 

d ' w 

partially  explained  by  investigating  these  two  derivatives. 
Both  and  are  functions  of  the  first  and  second 

derivatives  of  the  Prandtl  Number.  In  the  range  of 
temperatures  mentioned,  the  Prandtl  Number  as  a function 
of  temperature  reached  a maximum  value  and  also  contained 
an  inflection  point.  Therefore,  the  first  derivative  of 
the  Prandtl  Number  had  a sign  change,  and  both  derivatives 
varied  significantly  in  magnitude  in  this  range.  These 
changes,  coupled  with  changes  in  the  other  parameters  in 
the  p and  B derivatives,  caused  the  wall  temperature 
to  vary  as  shown  previously. 

Case  2 

This  case  involved  fixing  the  initial  conditions  and 
investigating  the  effects  of  the  fluctuations  of  the  five 
flow  parameters  independently.  The  initial  conditions  were 
fixed  as  the  baseline  conditions  given  previously  in  this 
chapter.  The  value  of  Qq  using  this  set  of  initial  condi- 
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tions  is  523.5  BTU/ft1 ' ^-sec . All  fluctuation  terms  were 
set  equal  to  zero  except  the  independent  variable  being 
fluctuated.  This  variable  was  fluctuated  over  a range  of 
0 < £ < 1 . 

For  the  entire  range,  fluctuations  of  the  density  at 
the  edge  of  the  boundary  layer  decreased  the  heat  transfer 
rate.  It  was  decreased  by  -0.3  percent  for  £.^=0.2  and  up 
to  -6.4  percent  for  £-^=1.0.  The  wall  density  fluctuations 
had  less  effect  on  the  heat  transfer  rate  than  any  of  the 
other  variables.  The  fluctuations  of  this  variable  decreased 
the  heat  transfer  rate  by  only  -4.5  percent  for  £2=1.0. 

Like  the  two  densities  previously  discussed,  the  stagnation 
pressure  fluctuations  also  decreased  the  heat  transfer  rate. 
It  was  decreased  by  -0.3  percent  for  £^=0.2  and  up  to  -9.4 
percent  for  £^=1.0. 

Fluctuations  in  the  temperature  at  the  wall  and  the 
temperature  at  the  edge  of  the  boundary  layer  had  a greater 
effect  on  the  heat  transfer  rate  than  did  the  other  three 
variables.  The  wall  temperature  fluctuations  decreased  the 
heat  transfer  rate  for  this  set  of  initial  conditions.  It 
was  decreased  by  -0.5  percent  for  a fluctuation  of  £^=0.2 
and  by  -12.8  percent  when  a fluctuation  of  £^=1.0  was 
considered.  Fluctuations  of  the  boundary  layer  temperature 
had  a greater  effect  on  the  heat  transfer  rate  than  did 
any  of  the  other  four  independent  variables.  Additionally, 
these  fluctuations  increased  the  heat  transfer  rate  in  con- 


trast to  the  other  four  variables.  The  fluctuations 


increased  the  heat  transfer  rate  by  0.?  percent  for  £.^=0.2 
and  up  to  17  • 5 percent  for  £~=1.0.  *he  results  of  the 
temperature  fluctuations  are  shown  in  Figures  10  and  11. 

The  results  of  three  other  specific  cases  are  shown  in 
Appendix  C . 
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Fig.  10.  Percent  Error  in  the  Heat  Transfer  Rate  for 
Fluctuations  of  Ts  (Baseline  Initial  Conditions) 
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Fig.  11.  Percent  Error  in  the  Heat  Transfer  Rate  for 
Fluctuations  of  Tw  (Baseline  Initial  Conditions) 


V . Conclusions 

The  effect  of  separate  and  simultaneous  fluctuations 
of  temperature  and  density  at  the  wall,  temperature  and 
density  at  the  edge  of  the  boundary  layer,  and  the  stagna- 
tion pressure,  on  the  heat  transfer  rate  as  calculated  by 
the  Fay-Riddell  equation,  have  been  theoretically  investi- 
gated. The  analysis  has  shown  that  neglecting  the  fluctua- 
tions of  these  variables  can  lead  to  substantial  errors  in 
the  calculated  heat  transfer  rate.  The  amount  of  the  error 
was  found  to  increase  with  the  magnitude  of  the  fluctuations. 

The  initial  value  of  the  temperature  at  the  edge  of 
the  boundary  layer  and  the  wall  temperature  had  a signifi- 
cant effect  on  the  amount  of  variation  of  the  heat  trans- 
fer rate.  The  percent  error  of  the  heat  transfer  rate  was 
decreased  as  the  initial  values  of  both  densities  and  the 
stagnation  pressure  were  increased  over  a wide  range  of 
initial  conditions,  but  the  magnitude  of  the  change  was 
very  small. 

The  trends  in  the  calculated  heat  transfer  rate,  as  a 
function  of  fluctuations  in  each  of  the  independent  variables, 
have  been  established.  It  has  been  shown  that  the  greatest 
error  was  produced  by  fluctuations  in  the  temperature  at 
the  edge  of  the  boundary  layer  when  all  five  variables 
were  uncorrelated  and  only  one  variable  fluctuated.  For 
the  baseline  set  of  initial  conditions  used  in  this 
investigation,  this  temperature  was  the  only  one  of  the 
five  variables  that  increased  the  heat  transfer  rate.  For 
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the  case  of  uncorrelation,  with  all  variables  allowed  to  5 

a 

fluctuate  simultaneously,  the  combined  negative  effect  of 

the  two  densities,  the  stagnation  pressure,  and  the  wall  ^ 

temperature  tended  to  dominate  the  positive  effect  of  the  j 

temperature  at  the  edge  of  the  boundary  layer  resulting  in 

a decrease  of  the  heat  transfer  rate.  Again,  the  amount  of 

decrease  was  dependent  upon  the  magnitude  of  the  fluctua-  j 

tions.  With  perfect  positive  correlation  and  all  variables 

allowed  to  fluctuate  simultaneously,  it  was  found  that  the  i 

j 

heat  transfer  rate  increased  for  all  cases  investigated. 

j 

It  was  pointed  out  that  certain  cross  terms  have  a much  \ 

greater  effect  on  the  heat  transfer  rate  than  others  ^ 

(i.e.  the  terms  involving  fluctuations  of  the  temperature  i 

1 

at  the  edge  of  the  boundary  layer).  j 

For  all  cases  studied,  the  trends  were  established  over 

j 

a wide  range  of  fluctuations  because  no  experimental  data 

was  available  on  the  magnitudes  of  the  fluctuations  or  1 

the  correlation  factors  for  the  cross  terms.  However,  ! 

when  a suitable  method  is  found  to  accurately  measure  or 
predict  the  fluctuations,  the  results  of  this  study  can 
be  applied  directly  to  show  the  effects  of  not  considering 
the  fluctuations  and  calculating  the  heat  transfer  rate 

i 

using  the  steady  laminar  Fay-Riddell  equation.  Additionally, 
the  significant  fluctuations  can  be  included  in  the  heat 
transfer  rate  calculations  by  the  same  modelling  technique 
used  in  this  study. 
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APPENDIX  A 

ANALYTICAL  DERIVATIVES  OF  THE  TRANSPORT 
AND  THERMODYNAMIC  PROPERTIES 
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APPENDIX  A 


Analytical  Derivatives  of  the  Transport 
and  Thermodynamic  Properties 


Enthalpy  Derivatives 

The  expression  used  for  the  species  enthalpy  was 


„ *rVT 

hj-hj  - 5+2(ry  |)  (r\,-\)  + £*i  9*.'  C 

T 2 TV'/l.i)  T|^e£vr 


( A-l ) 


Normalizing  the  temperature  with  respect  to  T ^ and  the 
energies  with  respect  to  RTre^  and  solving  for  the  species 
enthalpy 


- hj  + RT 
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where  the  temperature  and  energies  represent  normalized 
values.  The  first  derivative  of  the  enthalpy  with  respect 
to  temperature  was  expressed  as 
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The  second  derivative  of  the  enthalpy  with  respect  to 
temperature  was  obtained  by  taking  the  derivative  of  each 
term  in  the  brackets  seperately.  The  derivative  of  the 
first  term  was  zero.  The  derivative  of  the  second  was 
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expressed  as 
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After  simplifying  and  rearranging,  the  second  term  was 
written  as 
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The  derivative  of  the  third  term  in  Equation  A-3  involving 
the  electronic  levels  was 

T3  = [1 9,  [I  T 5,h)(^)(c'^)1 


(Wj/Tr 

Simplifying  the  first  part  of  this  expression  lead  to 
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The  derivative  of  the  fourth  term  in  Equation  A-3,  which 
also  involves  the  electronic  levels,  was 
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Simplifying  this  expression  lead  to 


t;  = 


Now  having  the  derivative  of  each  term,  the  second  deriva- 
tive of  the  enthalpy  was  expressed  as 

= WT,  (A- 10) 

Prandtl  Number  Derivatives 

The  Prandtl  Number  was  determined  using  curve  fit  data 
for  three  ranges  of  the  ratio  of  the  wall  enthalpy  and  a 
reference  enthalpy.  For  h ^/h  < 0.005.  the  Prandtl  Number 
was  set  equal  to  0.77.  Therefore,  for  this  range,  the 
first  and  second  derivatives  of  Pr  were  equal  to  zero. 

For  0 . 005  - hw/hQ £ 0 . 0 , the  Prandtl  Number  was  expressed 

by 
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where  P,  through  represented  the  coefficients  of 
Equation  20.  The  first  and  second  derivatives  of  the 
Prandtl  Number  were 
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For  0.1  £ h 2.0,  the  Prandtl  Number  was  expressed 


as 


Pr  = 


( A-14) 


where  P through  P1Q  represent  the  coefficients  of  Equation 
21.  The  first  and  second  derivatives  of  the  Prandtl 
Number  for  this  range  were 
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The  enthalpy  derivatives  were  shown  previously  in  this 
section . 


Viscosity  Derivatives 

The  viscosity  was  determined  by  the  Sutherland  Formula. 
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The  first  and  second  derivatives  of  the  viscosity  with 
respect  to  temperature  were 
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Lewis  Number  Derivatives 

The  Lewis  Number  was  determined  as  a function  of  the 

ratio  of  the  enthalpy  at  the  edge  of  the  boundary  layer  and 

a reference  enthalpy.  If  h /h  <0.1,  the  Lewis  Number  v/as 

s o 

set  equal  to  1.4.  Therefore,  the  first  and  second  deriva- 
tives of  the  Lewis  Number  were  set  equal  to  zero.  If 
hg/ho^  0.1,  the  Lewis  Number  was  expressed  by  Equation  22 
and  the  first  and  second  derivatives  were 
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h.  h/ 

The  enthalpy  derivatives  were  shown  previously  in  this 
section. 


(A-21) 
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Derivatives 

The  h^/hs  ratio  was  determined  using  curve  fit  data  as 
a function  of  the  enthalpy  at  the  edge  of  the  boundary  in 
BTU/lb.  The  ratio  was  expressed  as 

' (A-22) 

ns 

where  through  represented  the  coefficients  of  Equation 
23.  The  first  fid  second  derivatives  of  the  ratio  were 

(Hi-2  H3hs  + 3 H.+  hj  ' 4-  H5  Kg  ) / (A-23) 
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APPENDIX  B 


Computer  Program 
Explanation  of  Program 

The  Fortran  computer  program  incorporates  the  equations 
derived  in  Chapters  II  and  III  to  determine  the  effect  of 
fluctuations  of  the  five  independent  variables  on  the  heat 
transfer  rate  as  calculated  by  the  Fay-Riddell  equation. 

In  addition  to  the  main  program,  there  are  numerous 
subroutines . 

The  main  program  contains  the  calling  sequence  for  the 
various  subroutines.  The  heat  transfer  rate  without 
fluctuations  is  calculated  in  this  part  of  the  program. 

Also,  the  second  derivatives  that  do  not  involve  cross 
terms  are  computed  after  the  necessary  parameters  are 
determined  from  other  subroutines.  These  values  are 
printed  out  in  the  main  program. 

Subroutine  READ  initializes  the  five  independent 
variables  and  the  desired  fluctuations.  The  fluctuations 
are  normalized  by  the  initial  values  of  the  independent 
variables.  These  values  are  also  printed  out  in  this  sub- 
routine . 

Subroutine  THERM  computes  the  enthalpy  of  the  mixture 
after  computing  the  enthalpy  of  each  species.  Although 
only  the  nitrogen  and  oxygen  molecules  are  considered  in 
this  study,  the  subroutine  is  written  to  allow  for  easy 
inclusion  of  other  species.  The  temperatures  and  energies 
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are  normalized  for  the  enthalpy  computations.  The  enthalpy 
of  the  mixture  is  returned  to  the  main  program  in  the  units 
of  ft  /sec  . Subroutine  HDE  computes  the  first  and  second 
derivatives  of  the  enthalpy  in  the  same  manner.  The  temper- 
atures and  energies  are  also  normalized  to  calculate  these 
derivatives . 

Subroutine  VISCTY  determines  the  viscosity  at  the  wall 
and  the  boundary  layer  as  a function  of  the  respective 
temperatures.  The  viscosity  is  returned  to  the  main  pro- 
gram in  the  units  of  slugs/ft-sec . Subroutine  XMUDE  com- 
putes the  first  and  second  derivatives  of  the  viscosities. 

Subroutine  PRANUM  and  LEWIS  calculate  the  Prandtl 
Number  and  Lewis  Number  from  the  curve  fit  data  shown  in 
Chapter  III.  The  reference  enthalpy  used  in  both  sub- 
routines is  2.119x10®  ft^/sec^  (Ref  11).  Subroutines 
PRADE  and  XLEWDE  compute  the  first  and  second  derivatives 
of  the  Prandtl  Number  and  Lewis  Number,  respectively. 

Subroutine  RATIO  computes  the  ratio  of  the  dissociation 
energy  and  the  enthalpy  at  the  edge  of  the  boundary  layer. 
The  enthalpy  is  converted  to  BTU/lb  for  the  curve  fit  data 
used  in  this  subroutine.  The  first  and  second  derivatives 
of  this  ratio  are  computed  in  subroutine  HRATDE. 

Subroutines  ALPHDE , ADE,  BETDE,  and  BEEDE  calculate 
the  first  and  second  derivatives  of  » p A,  , and  B, 
respectively.  These  values  are  calculated  from  the  initial 
conditions  and  the  other  derivatives  previously  computed. 

The  magnitudes  of  the  derivatives  are  returned  to  the  main 


53 


program  and  used  to  compute  the  derivatives  of  the  heat  1 

transfer  rate  with  respect  to  the  temperature  at  the  wall  ^ 

and  the  temperature  at  the  edge  of  the  boundary  layer.  j 

Subroutines  GAWADE,  GASDE,  and  GAPTDE  compute  the 
first  and  second  derivatives  of  ^ with  respect  to  the 
density  at  the  wall,  density  at  the  edge  of  the  boundary 

layer,  and  the  stagnation  pressure.  The  values  of  these  j 

derivatives  are  returned  to  the  main  program  and  used  to 

compute  the  second  derivatives  of  the  heat  transfer  rate  j 

with  respect  to  the  same  three  independent  variables. 

Subroutines  CROSS  computes  all  of  the  second  derivatives 

j 

of  the  heat  transfer  rate  that  involve  cross  terms  of  the 

independent  variables.  The  value  of  these  derivatives  are  i 

| 

also  written  out  in  this  subroutine.  j 

J 

Subroutine  NORM  normalizes  and  non-dimensionalizes  all 

I 

j 

of  the  second  derivatives  of  the  heat  transfer  rate.  They 
are  normalized  by  dividing  each  term  by  the  magnitude  of  ' 

the  heat  transfer  rate  without  fluctuations.  The  values  of 
the  normalized  derivatives  are  then  written  out.  The  | 

derivatives  are  non-dimensionalized  by  multiplying  each 
second  derivative  by  the  magnitude  of  the  independent 
variable  or  variables  that  the  derivative  was  taken  with 
respect  to.  The  values  of  the  normalized,  non-dimensional- 
ized derivatives  are  also  printed  out  in  this  subroutine. 

In  subroutine  FLUCTS,  the  normalized,  non-dimensional- 
ized second  derivatives  of  the  heat  transfer  rate  are  multi- 
plied by  the  magnitude  of  the  desired  fluctuations.  These 
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terms  are  summed  to  obtain  the  heat  transfer  rate  from  the 


Taylor  Series  expansion.  The  exact  and  normalized  value  of 

« 

the  heat  transfer  rate  with  fluctuations  are  printed  out  in 
this  subroutine.  Additionally,  a check  is  conducted  to 
ensure  the  full  Taylor  Series  expansion  accurately 
approximates  the  exact  value  of  the  heat  transfer  rate  for 

the  fluctuations  considered.  1 

Fortran  Variables 

The  major  Fortran  variables  and  the  definition  of  each  \ 

are  listed  in  Table  II.  Additionally,  the  second  deriva-  I 

J 

tives  of  the  heat  transfer  rate  are  expressed  by  the  letters  j 

D2Q  followed  by  the  independent  variables  that  the  deriva-  j 

tive  was  taken  with  respect  to  (i.e.  the  second  derivative 

| 

of  the  heat  transfer  rate  with  respect  to  the  wall  tempera-  ! 

ture  and  density  is  D2QTWRW).  j 

The  first  and  second  derivatives  of  the  other  flow 
properties  have,  in  addition  to  the  letters  used  for  the 

flow  property  itself,  a D for  the  derivative  and  a 1 or  2 • 

I 

denoting  which  derivative  is  represented.  ; 

The  fluctuations  are  in  matrix  form  of  dimension  EPS 
(5,5).  The  diagonal  elements  represent  the  fluctuations 

of  only  one  variable  and  the  off-diagonal  elements  repre-  1 

sent  the  cross-term  fluctuations. 

j 

A listing  of  the  Fortran  program  follows  in  this 
Appendix. 
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Definitions  of  the  Major  Fortran  Variables 
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Variable  (Dimensions)  Definition  (Units) 


PROGRAM  H r AT  ( INPUT , TAPFS= I NPUT , OUTPUT, TA PE6=0UTPUT ) 
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RESULTS  OF  OTHER  SPECIFIC  CASES  INVESTIGATED 
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APPENDIX  C 


Results  of  Other  Specific  Cases  Investigated 

Case  3 

This  case  involved  fixing  the  initial  conditions  and 
investigating  the  effects  of  the  simultaneous  fluctuations 
of  the  five  independent  variables.  In  this  case  the 
fluctuations  are  assumed  to  be  uncorrelated.  The  five 
variables  are  fluctuated  simultaneously  over  a range  of 
Of  £ £ 1.  In  other  words  * 

= Ja/>  = (<£?>  = {a7>  = ^7> 

and 

<£,  a>  * o 

over  the  range  of  fluctuations  mentioned. 

The  results  of  this  investigation  are  shown  in  Figure 
12.  For  all  fluctuations  considered,  the  heat  transfer 
rate  decreased  and  the  magnitude  of  the  percent  error 
increased  as  the  value  of  the  fluctuations  increased.  The 
heat  transfer  rate  decreased  by  -0.6  percent  for  the  five 
independent  variables  fluctuated  for  £=0.2  and  decreased 
by  -15.8  percent  for  £ =1.0.  See  Figure  12  for  this  case. 

Case  4 

This  case  also  involved  fixing  the  initial  conditions 
and  investigating  the  effects  of  the  simultaneous  fluctua- 
tions of  the  five  independent  variables.  However,  in  this 
case,  perfect  positive  correlation  was  considered 
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PERCENT  ERROR  IN  Q 

-17. SO  -IS. 00  -12. SO  -10.00  -7.50  -5.00  -2. SO 


NORMAL  I ZED  FLUCTUATIONS 

§3.00  0.20  0.40  0.60  0.80 


1.00 


Q =523.5  BTU/ft1,5-sec 


Fig  12.  Percent  Error  in  the  Heat  Transfer  Rate  for 
Fluctuations  of  all  Variables  - Uncorrelated  (Base- 
line Initial  Conditions) 
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= \[<e7>  = ^7>  = 


={a7> 

and 

<elEiy  = ^7> 

The  range  of  fluctuations  was  again  0-£-  1.  The  re- 
sults of  this  investigation  are  shown  in  Figure  13 . The 
addition  of  the  cross  terms  signficantly  increased  the 
value  of  the  heat  transfer  rate.  The  heat  transfer  rate 
was  increased  by  2.5  percent  when  £ =0.2  and  was  increased 
up  to  63.^  percent  for  £=1.0.  This  was  a significant 
increase  in  the  heat  transfer  rate;  therefore,  the  effect 

of  the  cross  terms  is  investigated  further  in  the  next  case.  : 

Case  5 

This  case  involved  fixing  the  initial  conditions  and 
choosing  two  of  the  independent  variables  to  be  fluctuated. 

Positive  correlation  was  assumed  for  each  set  of  two 

# , * 
variables  chosen.  There  were  many  combinations  of  two 

variables  that  could  be  chosen.  The  magnitudes  of  the 

second  derivatives  of  the  heat  transfer  rate  used  in  the 

Taylor  Series  expansion  were  observed  to  determine  which 

terms  appeared  to  have  greater  effect  on  the  heat  transfer 

rate.  The  magnitudes  of  these  non-dimensionalized, 

normalized  derivatives  for  the  baseline  conditions  were 


*2°2  z 0-12  75 
eps 

d’°  = 0.03  75 

= -0.0900 

* lQ  - 0.  1 
<^TS 
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- 0.3^83 

A 5 

r -0.0097 

Ap„ATw 

= -O  2589 

^!q  =0.0250 

fp°a  = -0.187  5 

<9  lt2 

/r.Va03,T 

<7q 

s0;0  =0.0150 

«V°'3M5 

T^f-  = 0.2169 

A Ts 

= -0./4-6  3 

= -O.OZH-h- 

jTw}Pt, 

The  absolute  value  of  the  magnitudes  of  four  of  the  cross 
term  derivatives  was  greater  than  0.1.  These  derivatives 

Q anri  — o*  Q Since  three  of 


are 


^Ts«Jps  c^Ts^t* 

these  derivatives  involved  a fluctuation  in  the  temperature 

at  the  edge  of  the  boundary  layer,  all  four  cross  terms 
involving  this  variable  were  investigated.  Additionally, 
an  investigation  was  made  of  the  -- ^ term.  For  each 

& 6s  * *W 

case  checked,  the  two  variables  were  fluctuated  simultan- 
eously over  a range  of  0i  £ - 1. 


<£/£,>->/?£ ?>  /<£■/> 

All  of  the  cases  involving  fluctuations  of  the  temperature 
at  the  edge  of  the  boundary  layer  increased  the  heat  trans- 
fer rate.  When  it  was  fluctuated  with  the  density  at  the 
edge  of  the  boundary  layer,  the  heat  transfer  rate  was 
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increased  by  1.2  percent  for  2=0.2  and  up  to  32.7  percent 


for  8=1.0.  Similar  results  were  obtained  when  T was 

s 

fluctuated  with  the  wall  density.  The  results  of  these 
two  cases  are  shown  in  Figure  14. 

dw 

The  effect  on  the  heat  transfer  rate  was  the  largest 
when  the  temperature  at  the  edge  of  the  boundary  layer  was 
fluctuated  simultaneously  with  the  stagnation  pressure. 

The  heat  transfer  rate  was  increased  by  1.7  percent  for 
£=0.2  and  up  to  44.2  percent  for  £ = 1.0.  There  was  a 
much  smaller  effect  when  the  two  temperature  variables  were 
fluctuated  together.  The  results  of  these  two  cases  are 
shown  in  Figure  15 . 

Simultaneous  fluctuations  of  the  wall  temperature  and 
the  density  at  the  edge  of  the  boundary  layer  revealed  a 
significant  decrease  in  the  heat  transfer  rate.  It  was 
decreased  by  -0.9  percent  for  £=0.2  and  by  -20.8  percent 
for  £=1.0.  Results  of  this  case  are  shown  in  Figure  16. 

In  addition  to  investigating  these  five  cases,  some 
interesting  results  were  obtained  from  investigating  the 
analytical  expressions  of  the  various  second  derivatives  in 
the  Taylor  Series  expansion.  Inspection  of  the  cross  terms 
involving  combinations  of  the  two  densities  and  the  stagna- 
tion pressure  revealed  that  the  magnitude  of  these  second 
derivatives  could  be  easily  obtained.  The  three  second 
derivatives  are  simply  the  product  of  a constant  and  the 
heat  transfer  rate  without  fluctuations.  Although  those 
terms  have  a small  effect  compared  to  some  of  the  other 
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cross  terms,  inclusion  of  them  into  the  Taylor  Series 
expansion  is  easily  accomplished  if  desired.  An  inspection 
of  any  of  the  analytical  derivatives  in  the  Taylor  Series 
expansion  can  lead  to  conclusions  about  their  dependency 
upon  the  transport  and  thermodynamic  properties  as  well 
as  the  other  flow  parameters. 
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PERCENT  ERROR  IN 
cP.OO  5.00  10.00  15.00  20.00 
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□ Fluctuations  of  Ts  and  q 
O Fluctuations  of  Ts  and  p 
Q =523.5  BTU/ft1* 5-sec 


NORMALIZED  FLUCTUATIONS 


Fig.  14.  Percent  Error  in  the  Heat  Transfer  Rate  for 
Fluctuations  of  (a)  Ts  and  <?s,  and  (b)  Ts  and  pw  - 
Correlated  (Baseline  Initial  Conditions) 


o 
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NORMALIZED  FLUCTUATIONS 


Fig.  15.  Percent  Error  in  the  Heat  Transfer  Rate  for 
Fluctuations  of  (a)  Ts  and  ^to.  and  (b)  Ts  and  Tw  - 
Correlated  (Baseline  Initial  Conditions) 
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